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Some cylindrically symmetric inhomogeneous viscous fluid cosmological models with electro- 
magnetic field are obtained. To get a solution a supplementary condition between metric potentials 
is used. The viscosity coefficient of bulk viscous fluid is assumed to be a power function of mass 
density. Without assuming any ad hoc law, we obtain a cosmological constant as a decreasing function 
of time. The behaviour of the electro-magnetic field tensor together with some physical aspects of the 
model are also discussed. 
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1. Introduction 

Inhomogeneous cosmological models play an important role in understanding some es- 
sential features of the universe such as the formation of galaxies during the early stages of 
evolution and process of homogenization. The early attempts at the construction of such 
models have done by Tolman 1 and Bondi 2 who considered spherically symmetric mod- 
els. Inhomogeneous plane-symmetric models were considered by Taub 3,4 and later by 
Tomimura, 5 Szekeres, 6 Collins and Szafron, 7 Szafron and Collins. 8 Recently, Senovilla 9 
obtained a new class of exact solutions of Einstein's equation without big bang singu- 
larity, representing a cylindrically symmetric, inhomogeneous cosmological model filled 
with perfect fluid which is smooth and regular everywhere satisfying energy and causality 
conditions. Later, Ruis and Senovilla 10 have separated out a fairly large class of singular- 
ity free models through a comprehensive study of general cylindrically symmetric metric 
with separable function of r and t as metric coefficients. Dadhich et al. 11 have estab- 
lished a link between the FRW model and the singularity free family by deducing the 
latter through a natural and simple inhomogenization and anisotropization of the former. 
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Recently, Patel et al. 12 presented a general class of inhomogeneous cosmological models 
filled with non-thermalized perfect fluid by assuming that the background spacetime ad- 
mits two space-likecommuting killing vectors and has separable metric coefficients. Bali 
and Tyagi 13 obtained a plane-symmetric inhomogeneous cosmological models of perfect 
fluid distribution with electro-magnetic field. Recently, Pradhan et al. 14 have investigated a 
plane-symmetric inhomogeneous viscous fluid cosmological models with electro-magnetic 
field. 

Models with a relic cosmological constant A have received considerable attention re- 
cently among researchers for various reasons (see Refs. 15 ~ 19 and references therein). Some 
of the recent discussions on the cosmological constant "problem" and consequence on cos- 
mology with a time-varying cosmological constant by Ratra and Peebles, 20 Dolgov 21-23 
and Sahni and Starobinsky 24 have pointed out that in the absence of any interaction with 
matter or radiation, the cosmological constant remains a "constant". However, in the pres- 
ence of interactions with matter or radiation, a solution of Einstein equations and the as- 
sumed equation of covariant conservation of stress-energy with a time-varying A can be 
found. For these solutions, conservation of energy requires decrease in the energy den- 
sity of the vacuum component to be compensated by a corresponding increase in the en- 
ergy density of matter or radiation. Earlier researchers on this topic, are contained in 
Zeldovich, 25 Weinberg 16 and Carroll, Press and Turner. 26 Recent observations by Perl- 
mutter et al. 27 and Riess et al. 28 strongly favour a significant and positive value of A. 
Their finding arise from the study of more than 50 type la supernovae with redshifts in 
the range 0.10 < z < 0.83 and these suggest Friedmann models with negative pressure 
matter such as a cosmological constant (A), domain walls or cosmic strings (Vilenkin, 29 
Garnavich et al. 30 ) Recently, Carmeli and Kuzmenko 31 have shown that the cosmologi- 
cal relativistic theory (Behar and Carmeli 32 ) predicts the value for cosmological constant 
A = 1.934 x 10~ 35 s~ 2 . This value of "A" is in excellent agreement with the measurements 
recently obtained by the High-Z Supernova Team and Supernova Cosmological Project 
(Garnavich et al. 30 ; Perlmutter et al. 27 ; Riess et al. 28 ; Schmidt et al. 31 ) The main conclu- 
sion of these observations is that the expansion of the universe is accelerating. 

Several ansatz have been proposed in which the A term decays with time (see Refs. 
Gasperini, 34,35 Berman, 36 Freese et al., 19 Ozer and Taha, 19 Peebles and Ratra, 37 Chen 
and Hu, 38 Abdussattar and Viswakarma, 39 Gariel and Le Denmat, 40 Pradhan et al. 41 ). Of 
the special interest is the ansatz A oc S~ 2 (where S is the scale factor of the Robertson- 
Walker metric) by Chen and Wu, 38 which has been considered/modified by several authors 
( Abdel-Rahaman, 42 Carvalho et al., 19 Waga, 43 Silveira and Waga, 44 Vishwakarma 45 ). 

Most cosmological models assume that the matter in the universe can be described 
by 'dust' (a pressure-less distribution) or at best a perfect fluid. However, bulk viscosity 
is expected to play an important role at certain stages of expanding universe. 46-48 It has 
been shown that bulk viscosity leads to inflationary like solution, 49 and acts like a negative 
energy field in an expanding universe. 50 Furthermore, there are several processes which 
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are expected to give rise to viscous effects. These are the decoupling of neutrinos during 
the radiation era and the decoupling of radiation and matter during the recombination era. 
Bulk viscosity is associated with the Grand Unification Theories (GUT) phase transition 
and string creation. Thus, we should consider the presence of a material distribution other 
than a perfect fluid to have realistic cosmological models (see Gr0n 51 for a review on cos- 
mological models with bulk viscosity). A number of authors have discussed cosmological 
solutions with bulk viscosity in various context. 51-54 

The occurrence of magnetic fields on galactic scale is well-established fact today, and 
their importance for a variety of astrophysical phenomena is generally acknowledged as 
pointed out Zeldovich et al. 55 Also Harrison 56 has suggested that magnetic field could 
have a cosmological origin. As a natural consequences, we should include magnetic fields 
in the energy-momentum tensor of the early universe. The choice of anisotropic cosmolog- 
ical models in Einstein system of field equations leads to the cosmological models more 
general than Robertson-Walker model. 57 The presence of primordial magnetic fields in the 
early stages of the evolution of the universe has been discussed by several authors. 58 ~ 67 
Strong magnetic fields can be created due to adiabatic compression in clusters of galax- 
ies. Large-scale magnetic fields give rise to anisotropies in the universe. The anisotropic 
pressure created by the magnetic fields dominates the evolution of the shear anisotropy and 
it decays slower than if the pressure was isotropic. 68 ' 69 Such fields can be generated at 
the end of an inflationary epoch. 70 ~ 74 Anisotropic magnetic field models have significant 
contribution in the evolution of galaxies and stellar objects. Bali and Ali 75 had obtained 
a magnetized cylindrically symmetric universe with an electrically neutral perfect fluid as 
the source of matter. Several authors 76-81 have investigated Bianchi type I cosmological 
models with a magnetic field in different context. 

Recently Singh et al. 82 obtained some cylindrically symmetric inhomogeneous cosmo- 
logical models for a perfect fluid distribution with electro-magnetic field. Motivated the 
situations discussed above, in this paper, we shall focus upon the problem of establishing a 
formalism for studying the general relativistic evolution magnetic inhomogeneities in pres- 
ence of bulk viscous in an expanding universe. We do this by extending the work of Singh 
et al. 82 by including an electrically neutral bulk viscous fluid as the source of matter in 
the energy-momentum tensor. This paper is organised as follows. The metric and the field 
equations are presented in section 2. In section 3 we deal with the solution of the field 
equations in presence of bulk viscous fluid. The sections 3.1, 3.2 and 3.3 contain the three 
cases as K > 0, K < and K = and also contain some physical aspects of these models 
respectively. Finally in section 4 concluding remarks have been given. 



2. The metric and field equations 

We consider the metric in the form 

ds 2 = A 2 (dx 2 - dt 2 ) + B 2 dy 2 + C 2 dz 2 , 



(1) 
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where the metric potential A is a function of t alone and B and C are functions of x and t. 
The energy momentum tensor in the presence of bulk stress has the form 



TV = {p + p) Vi v j + pgf+Ef, 
where Ef is the electro-magnetic field given by Lichnerowicz 83 as 

1 



E 



p 



\h\ 2 Vi v> 



KM 



and 



(2) 



(3) 



(4) 



p = p- &.i 

Here p, p, p and £ are the energy density, isotropic pressure , effective pressure, bulk viscous 
coefficient respectively and v % is the flow vector satisfying the relation 



gijv'v 1 = -1 

p is the magnetic permeability and hi the magnetic flux vector defined by 

hi = 4 *F iiV j 

where *Fij is the dual electro-magnetic field tensor defined by Synge 84 to be 

V=9 



" p. 



-tijkl-F 



(5) 



(6) 



(7) 



Fij is the electro-magnetic field tensor and is the Levi-Civita tensor density. Here, 
the comoving coordinates are taken to be v 1 = = v 2 = v 3 and v 4 = -j. We consider the 
current to be flowing along the z-axis so that F12 is the only non- vanishing component of 
F^. The Maxwell's equations 

= 0, (8) 



J 1 



(9) 



require that F i2 is a function of x alone. Here the semicolon represents a covariant differ- 
entiation. We assume that the magnetic permeability is a function of x and t. The Einstein's 
field equations ( in gravitational units c = 1, G = 1 ) read as 



(10) 



for the line element (1) has been set up as 

87tA 2 f p - 



F? 2 



A4B4 A4C4 



AB 



AC 
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(11) 
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where 



0„ A2 ( - i F 12 \ A 4 ^44 C44 Cll . ,2 n ~, 

87TA { P+ 2^)=^- — - — + — - AA > (12) 

Z^A 2 f* F l 2 2 \ A 4 A 44 B44 B n . , 2 , nl 



2flA 2 B 2 

,t« .ir b c A ' 

B14 Ci, 





B4C4 

1 — 








d\ 





Sl4 + d4_^ (15) 



£ , B 4 C 4 



The suffixes 1 and 4 by the symbols A, B and C denote differentiation with respect to x 
and t respectively. 



3. Solution of the field equations 

From Equations (1 1), (12) and (13), we have 

A44 A\ A4B4 A4C4 B44 B4C4 Cu B\C\ _ 
~A~~A^ ~AB ~AC B BC (T ~BC~ ~ { ' 

8ttFi2 _ C44 Cu B44 Bn 

p,B 2 ~ C C B B ( ' 

Equations (1 1) - (15) represent a system of five equations in eight unknowns A, B, C, p, p, 
F12, A and £. To get a determinate solution, we need three extra conditions. 
Firstly following Singh et a/. 81 , we assume a supplementary condition between the metric 
potentials as 

A4 ( B4 C4 



(18) 



and 

Cn _ B1C1 

C BC ' ' 
where n and K are arbitrary constants. Let us consider that 

B = f(x)g(t) 

C = h{x)k{t) (20) 
Using (18) and (20) in Equations (15), (16) and (17), we have 
(n + l)i_ n £ /y/ 



(l-n)£ + n£ h'/h 



= a (constant, say), (21) 
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{n - 1) - 9 - n k-Tk = ' (22) 

where prime and dot stand for differentiation with respect to x and t respectively. Integrat- 
ing (21) leads to 

/ = Lh- a , 

g = M 1 k^^ 1 (24) 

where L and Mi are constants of integration. 

Here we consider three cases according to the values of K. 



3.1. Case (1): K > 

Using (24) in Equations (22) and (23) and after making suitable transformation of co- 
ordinates, the geometry of the spacetime (1) reduces to the form 

ds 2 = a 2 cosher 11 (K^idX 2 - dT 2 ) 

+ cosh-^TT(K 2 X) cosh^^^ (K^dY 2 
+ cos\^{K 2 X)cosYr^ ±1 (K l T)dZ 2 , (25) 

where 

Kx = VKN, 
K 2 = y/K(a+l), 
N = 2na -2n-a-l, 

n(l-o) 

a = LMxa^"^" . 
The effective pressure and density for the model (25) are given by 

8np = 8tt(p - £0) = 



a 2 cosh ™ (KiT) 

(, °? K T K 2 t&nh 2 (K 1 T)\ ( 

(n - na - a + 2)K + tanh 2 (K 2 ^) + K 3 v ; - A, 

\ Qt \ 1 1 / 



(26) 



87r P = 2n(c-l) 

a 2 cosh « (KiT) 



. a(a-2)K l2 ,„^ T ^ K 2 tanh' ! (K 1 T)\ A 
{na-n- a)K - y — ^ t&nh 2 (K 2 X) + K 3 ^ ' j + A, (27) 
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where 

K 3 = n 2 {3a 2 - Aa + 2) - 2n(a - l)(a - 3) - 2(a 2 - a - 1). 



Here 6* is the scalar of expansion calculated for the flow vector v % as 

AT^tanh^iT) 



iViif tanh(ifiT) 



where 

iVi = 3na — 3n — a — 1 
For the specification of £, we assume that the fluid obeys an equation of state of the form 

P = IP, (29) 

where 7(0 < 7 < 1) is a constant. 

Thus, given we can solve for the cosmological parameters. In most of the investi- 
gations involving bulk viscosity is assumed to be a simple power function of the energy 
density. 85 " 87 

£(*) = SaP m , (30) 

where £ an d m are constants. If m = 1, Eq. (30) may correspond to a radiative fluid. 88 
However, more realistic models 89 are based on m lying in the regime < m < \. 
On using (30) in (26), we obtain 

8tt(p - ^p m 6) = ^ x 

a 2 cosh » {K\T) 

((»-««-« + 2)K + ^ t,nh 2 (K 2 X) + K ^t^ 2 { Kl T) ^ ^ ^ 

3.1.1. Model I: (£ = £0) 

When m = 0, Equation (30) reduces to £ = £ - With the use of Equations (27), (28) 
and (29), Equation (31) reduces to 

8tt(1 + 7)p = + 

aifi cosh - ^ (KiT) 

2(1 - a) (a + \)KK\ + 2aK K 2 t&nh 2 (K 2 X) + 2(a + l)if 2 if 3 tanh 2 (XiT) 
(a + l)a 2 X 2 cosh^^ 11 (XiT) 
Eliminating p(t) between (27) and (32), we get 

(1 + 7)A = ^rry h 

aATi cosh^^ - {K{T) 
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X 



[{n — na — a + 2 — (na — n — a)^}(a + 1)KK 1 + 
{a + (a - 2)i\aKK\ tanh 2 (K 2 ^) 
+ (1 - 7) (a + l)K 2 K 2 tanh 2 (XiT)] 
3.1.2. Model II: (£ = £ p) 



(33) 



When m = 1, Equation (30) reduces to £ = £ P- With the use of (27), (28) and (29), 
Equation (31) reduces to 



8irp 



1+7 



2(1 - c?)KK\ + laKKl tanh 2 (iT 2 X) + 2 (a + l)K 2 K 3 tanh 2 (^iT) 



2n(a-l) 



(a + l)a 2 i^ cosh « (XiT) 
Eliminating between (27) and (34), we get 



(34) 



A 



1 + 7- 



N 1 ^ Kta,iih(K 1 T) 
aKi cosh 111 ^ [K\T) 



(n-na-a + 2) (a + l)ifi\: 2 + a 2 KK\ tanh 2 (ivT 2 ^) + (a + l)if 2 if 3 tanh^i^T) 

(a + l)a 2 if 2 cosh 211 ^ 11 (K\T) 
(na -n - a) (a + 1)KK 2 - a(a - 2)KK 2 t&nh 2 (K 2 X) + (a + l)K 2 K 3 tanh 2 (KiT) 

2n(o-l) 

,2 ^2 „^ c v, 5v 



(a + l)a 2 if 2 cosh^— (KiT) 
N 1 £, Kt<mh(K 1 T) 



7 



(35) 



From Equations (33) and (35), we observe that the cosmological constant is a decreasing 
function of time and it approaches a small value as time progresses (i.e., the present epoch), 
which explains the small value of A at present. Figure 1 clearly shows this behaviour of A 
as decreasing function of time in both the models I and II. 



Some Physical Aspects of the Models 

We shall now give the expressions for kinematical quantities and the components of con- 
formal curvature tensor. With regard to the kinematical properties of the velocity vector v l 
in the metric (25), a straight forward calculation leads to the following expressions for the 
shear of the fluid. 



ffn = a(Q 3 ^ i 1) ^ cosh I ^(X 1 r)tanh(^ 1 T), 



(36) 
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Figure 1 : The plot of cosmological constant (A) with time for model K > with param- 
eters a — 1.5, 7 = 0.5, K = 1, a = 1, £o = 1. solid line for power index m = 0, and 
dashed line for m — 1 respectively. 

0-22 = ^^r) K cosher (K 2 X)cosh : ' S ^ 1 (K 1 T)tan]i(K 1 T), (37) 
6aK\ 

^33 = - n 2 "}^ cosher (K 2 X) cosri"°~«~ 2 ° (i^T) tanh(^iT), (38) 

cr 44 = ^— L_ cosh 11 ^ (i^T) tanh(^iT), (39) 

where all other components, the rotation tensor lo and acceleration vanish. The expansion 
scalar 9 has already been given by (49). Since Vi = Vi y-u 3 = 0, the motion is geodetic. 
The non-vanishing physical components of conformal curvature tensor are given by 

C\l = CH = ^-^ [AaKt&nh 2 {K 2 X) - (a + 1)K+ 

6a 2 cosh^^(isfiT) 
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{n 2 (a - l) 2 - n(5a 2 - 4a + 1) + 2a} 
1 



if 2 tanh 2 (if iT) 



r <13 _ /-»24 
°13 — °24 



6a 2 cosh'"™ (KiT) 



[{K-2aKla,xh 2 (K 2 X) + 



(40) 



/^14 _ ,-y23 _ 
°14 — °23 — 



{2n(a-2)(a-l) + 2a} 
1 



if 2 tanh 2 (i4riT) 



6a 2 cosh^ 11 ^ 1 {KiT) 



2{na 2 -2a- l)K 2 



(41) 



/-ill /~i 

°12 — °13 



tanh^(^iT)- 

L 

2a.fi: tanh 2 (K 2 X) - 2K] (42) 
., _ n(l - a)V(« + l)^ tanh^T) tanh(M) 



2a 2 i^i cosh 22 ^^^) 



The non-vanishing component F12 of electro-magnetic field tensor and J 2 , the component 
of charge current density, are obtained as 



F? 2 



UaK\ 2(2o + l) JV + a-l 

r \ sech^T^(^ 2 X)cosh^^(if 1 T) 



47T 



(44) 



1 



4a 2 V ^7r(a + 1) 



cosh 



\K 2 X) sech~^ {K{T) 



p!^f{a + y) - 2a^fita,nh(K 2 X) 



(45) 



The source of electro-magnetic field exist, when matter and charge are present and does not 
exist when matter is absent. The models represent expanding, shearing, non-rotating and 
Petrov type I non-degenerate in general, in which the flow is geodetic. The model starts 
expanding at T > and goes on expanding indefinitely. It is observed that the expansion is 
minimum at T = 0. Since lim T ^oo f 7^ 0, the models do not approach isotropy for large 
values of T. The pressure contrast ^ and density contrast — tends to zero for large values 
of T which shows that inhomogeneity dies out for large values of T. When a = —1 the 
model reduces to its homogeneous form. It is remarkable to mention here that at the time 
of minimum expansion matter density dominates over the expansion in the model whereas 
at the time of maximum expansion, expansion in the model dominates over matter density. 



3.2. Case (2) : K < 0, let K = -M, M > 

Using (24) in Equations (22) and (23) and after making suitable transformation of co- 
ordinates, the geometry of the spacetime (1) reduces to the form 

ds 2 = b 2 cos^r 11 (K A T){dX 2 - dT 2 ) 

+ cos-^(K 5 X) cos- (W+ r 1)) (K 4 T)dY 2 
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2 -(N ++CX — 1) 

+ cos— (K 5 X) cos * (K 4 T)dZ 2 , 



where 



(46) 



K 4 = VMN, 
K 5 = y/M(a+l), 

n(l- a ) 

b = LM 1 a' 2 n - na +" 
The effective pressure and density for the model (46) are given by 

8-Kp = 8tt( p - £0) = ' 



b 2 cos^r 21 (K 4 T) ' 



/»r 2 " 2M 2,t, ^ ^ M 2 tan 2 (K 4 T) 

-(N+1)M+ - ; - tan 2 (if 5 X) + jf 6 - 



a + 1 

87T/3 = 



A, 



1 



6 2 cos 2l ^(K4T) 



(2-a)M- 



a(a-2)M 4 _ 2(r ^ ^ , ^ M 2 tan 2 (if 4 T) 



a + 1 



tan 2 (X 5 X) +if 7 - 



+ A, 



where 



K% = n(a — l)(5n — 5na + 4a + 4) — a 2 — a — 1, 
if 7 = n(a - l)(3na - 3n - 4a) + a 2 + a - 1. 
and 6 is the scalar of expansion calculated for the flow vector v l as 

d _ A^MtanfT^T) 



(47) 



(48) 



(49) 



3.24. Model I: (£ = £o) 

When m = 0, Equation (30) reduces to £ = £ - in this case Equation (47) with the use 
of (29), (48) and (49) reduces to 



8tt(1 +7)p : 



bKiCos 11 ^ {K 4 T) 



-(N + l)(a + 1)M Jf| + a(a + 2)MK 2 t&n 2 (K 5 X) + (K e + K 7 ){a + l)M 2 t&n 2 {K 4 T) 



(a + l)b 2 K 2 cos^tt^ 1 (K 4 T) 
Eliminating p(t) between (48) and (50), we get 



(50) 



(1 + 7)A = 



bK 4 costT 1 ^ 1 (K 4 T) 
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— X 



2n(l-o) 



(a + l)b 2 Kl cos "" « (if 4 T) 
[{2(n - 7) - a(2n - 7) + a}(a + l)MK_ 
{a 2 (2 + 7) - 2a-/} MKj tan 2 (K 5 X) 
+{K 6 - jK 7 )(a + l)M 2 t&n 2 {K 4 T)] 



(51) 



3.2.2. Model II: (£ = £op) 

When m = 1, Equation (30) reduces to £ 
use of (29), (48) and (49), reduces to 



top. In this case Equation (47) with the 



87T/9 



1+7- 



N^M t&n{K 4 T) 
bK 4 coa*^ (K A T)_ 



2(n+ l)(a 2 - l)MK 2 + a(a + 2)MKjttm 2 (K 5 X) + (K 6 + K 7 ){a + 1)M 2 tan 2 (if 4 T) 



(a + l)b 2 K 2 cos 
Eliminating p(t) between (48) and (52), we get 



^{K 4 T) 



(52) 



1+7- 



N lc :oMta,Ti(K 4 T) 
bK 4 cos :li ^ 1 {K 4 T)_ 



(2n - 2na - a) (a + \)MK\ + 2u 2 MK 2 4 tan 2 (K 5 X) + (a + l)K 6 M 2 t&n 2 (K 4 T) 



(a+ l)& 2 iff cos' 



L (if 4 T) 



(2 - a){a + \)MK\ - a(a - 2)MK\ tan 2 (K 5 X) + (a + 1)M 2 K 7 tan 2 (K 4 T) 



(a + l)b 2 K 2 cos' 



' (K4T) 



7' 



N&M tan(K 4 T) 
bK 4 cos 1 ^ 1 (K 4 T) 



(53) 



Our study for the case K < shows constant value of cosmological constant (A) for large 
values of time and do not decrease with time (this means that the universe in not expanding 
or may be steady state condition). In this case, detailed study shows that the scalar of ex- 
pansion 6 does not increase with time. Our study is inconsistent with work done by Singh 
et al. 82 The Singh et al. claim that the universe is expanding which does not match with 
our result. Also the claim of minimum and maximum expansion rate in 9 is reflection of 
periodicity of trigonometric functions involved there. We are trying to find feasible inter- 
pretation and situations relevant to K < 0. Further study is in progress. 

Some Physical Aspects of the Models 

The non-vanishing components of shear tensor cry are obtained as 



(Til = 



b(a 



3K 4 



1)M 

sec 



■{K 4 T) tan(if 4 T), 



(54) 
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J -2-> = ^ OL r) M SCC ^ i K 5 x ) sec 2 ^"^ {K 4 T) tan(K 4 T), (55) 



CT33 



3bK 4 
(1 - 2a)M 

ZbK 4 



sec~ — (K 5 X) sec "" ~ (.K4T) tan(i<: 4 T) , (56) 



CT 44 - sec^^GfC^tan^T) 
3 0X4 



(57) 



The expansion scalar 9 has already been given by (49). The rotation u is identically zero. 
Since tij — Vi -jV^ = 0, the motion is geodetic. 

The non-vanishing physical components of conformal curvature tensor are given by 



^12 _ ^34 

°12 — °34 



1 



66 2 cos 22 ^ (^T) 



[(4a + n - na)M + 4aM tan 2 (i^ 5 X) + 



{n(a - l)(3n - 3na + 5a - 1) - 4a 2 + 2} 



M tan 2 (K4T) 



N 



(58) 



r»13 _ ^24 _ 
°13 — °24 — 



66 2 cos^^tt^ (K 4 T) 



{2n(l-a)(a-2)+a 2 -2}__ 2 . 



*2 



-2aMsec 2 (K 5 X)] 



M 2 tan 2 (if 4 T) 
(59) 



-ol4 _ n 23 
o 14 — (_/< 



14 — ^23 — " 2»(i-«) , " 

6b 2 cos^r- 1 (K 4 T) 



2(a 2 — na 2 + n + 1) , r9 , . r ^ m . 

^^M 2 tan 2 (X 4 T) 



-2aM scc 2 (K 5 X)] 



, 2I _ n(a - l)y/(a + 1)M* tan(if 4 T) tan(jf 5 ^Q 

2& 2 ^ 4 cos 2 " ( » a) (K 4 T) 



r<2A _ s~t34 _ 
°12 — °13 — 



(60) 
(61) 



The non- vanishing component F 12 of electro-magnetic field tensor and J 2 , the component 
of charge current density, are obtained as 

2 = ^ 



12 



\ 47T J 

{a 2 (2n - 1) - 2n(2a - 1) + 1} 



sec^+r (K 5 X)sec^^- {K 4 T) 

Mt&n 2 (K 4 T) -at<m 2 (K 5 X) - 1 



(62) 



M 



-<"i+°) , 
sech — " — (K 4 T) 



where 



#1 



46 2 y finHtfa + l) 
sec^TT{K 5 X) hfia{a + 1)Vm taa(K 5 X) 
sec 2 (K 5 X) - p,'H iy /(a + 1) + 2/zifi VM t&n{K 5 X) 

{a 2 (2n-l)-2n(2a-l) + ll . 2 



(63) 



M 2 tan 2 (K 4 T) - atan 2 (K 5 X) - 1 
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The models represent shearing, non-rotating and Petrol type I non-degenerate in general, in 
which the flow is geodetic. The model starts expanding at T > but the initial expansion 
is slow. When T is closer to -^q, it has stiff rise in the expansion then decreases. This 
shows the case of oscillation. It is observed that the expansion is minimum at T = or 
T = The large values of 9 near T — is reflection of trigonometric property. 
But expansion remains finite. Since liniT— oo § ^ 0, hence the models do not approach 
isotropy for large values of T. In this case inhomogeneity also dies out for large value of T. 

3.3. Case (3): K = 

Using (24) in Equations (22) and (23) and after making suitable transformation of co- 
ordinates, the geometry of the spacetime (1) reduces to the form 

2n(a — 1) , ~ -2a N + a-1 2 N-a + 1 n 

ds 2 = c 2 T^—{dX 2 - dT 2 )+X—T—^dY 2 +X — T^—dZ 2 , (64) 
where 

.. n(l-q) 

c = LM1CI2 . 

It is observed that the model (64) starts to expand from its singularity stage i.e. at T = 
and goes expanding indefinitely when T — > 00. The model represents expanding, shearing, 
non-rotating and Petrov type I non- degenerate in general in which the flow is geodetic. 
This model also does not approach isotropy for large values of T. It is also observed here 
at the time of initial singularity the matter density dominates over the expansion of the 
model. 

4. Conclusions 

We have obtained a new class of cylindrically symmetric inhomogeneous cosmological 
models electro-magnetic bulk viscous fluid as the source of matter. Generally the models 
represent expanding, shearing, non-rotating and Petrov type-I non-degenerate in which the 
flow vector is geodetic. In all these models, we observe that they do not approach isotropy 
for large values of time. It is concluded that (i) if a > then for the cases K > or the 
values of K tending to 0, the models have point type singularity at the time of maximum 
expansion and (ii) if < a < 1 then they have infinite singularity. Whereas the case 
K < is just opposite of the case K > 0. In all the cases the spacetime is conformally flat 
for large values of T and at the time of minimum expansion matter density dominates over 
expansion. For K = the material energy is more dominant over magnetic energy. 
The cosmological constant in all models given in section 3 are decreasing function of time 
and they all approach a small value as time increases (i.e., the present epoch) except the 
case K < . The values of cosmological "constant" for these models are found to be 
small and positive which are supported by the results from recent supernovae observations 
recently obtained by the High - z Supernova Team and Supernova Cosmological Project ( 
Garnavich et al. 30 ; Perlmutter et al. 27 ; Riess et a/. 28 ; Schmidt et al. 33 ). Thus, with our 
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approach, we obtain a physically relevant decay law for the cosmological constant unlike 
other investigators where adhoc laws were used to arrive at a mathematical expressions for 
the decaying vacuum energy. Thus our models are more general than those studied earlier. 
We would like to find feasible situations for K < 0. Our strong point of this model is 
that it in-cooperates matter density naturally and so makes feasible model which can be 
in-cooperates the physical constraints. 
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